program vectortest
	implicit none
	double precision:: a(3), b
	a(1) = 1d0
	a(2) = 2d0
	a(3) = 3d0
	b = sum( a**2 )
	print *, b
end program
